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Abstract 

The theory of NQR spin diffusion is extended to the case of spin lattice relaxation and spin 
diffusion in an inhomogeneous field. Two coupled equations describing the mutual relaxation and 
the spin diffusion of the nuclear magnetization and dipolar energy were obtained by using the 
method of nonequilibrium state operator. The equations were solved for short and long times 
approximation corresponding to the direct and diffusion relaxation regimes. 



1 



Introduction 

Studies of tlie NMR |l[ S Ql and NQR of nuclei have demonstrated that spin diffusion 
plays an important role in the relaxation of nuclei in the presence of paramagnetic impurities 
(PI). Such type of relaxation originates from the magnetic dipole-dipole interaction of PI 
with neighboring nuclei, which is inversely proportional to the sixth power of the distance. 
Thus, near the PI the equilibrium with the lattice is reached at a faster rate 
The nuclear magnetization during the relaxation process is spatially inhomogeneous over 
a sample volume, and this induces a spatial diffusion of the nuclear spin energy by, for 
example, flip-flop transitions due to the dipole-dipole interactions between nuclear spin. 

However, first, till now most studies of the nuclear spin diffusion were related to systems 
with nuclear spin I = 1/2, described by the Hamiltonian, whose main term includes just 
linear functions of spin operators and, correspondingly, forms equidistant energy spectra and, 
second, most of them deal with the process of a spin diffusion in homogeneous magnetic fields 

In many samples spin systems consist of nuclei with / > 1/2 and they interact with their 
environment through the electric quadrupole moment Q, and these interactions are strong 
enough to observe magnetic resonance of nuclei in the absence of an external magnetic field 
(pure NQR-case). Unlike the NMR-case, the NQR energy spectrum is non-equidistant and, 
in many cases, degenerated. These circumstances lead to certain difficulties in obtaining a 
diffusion equation and a calculation of a diffusion coefficient. 

Bisides, the spin diffusion processes no longer exactly conserve nuclear quadrupole energy 
in an inhomogeneous field because the quadrupole interaction energy is not identical for 
neighboring nuclear spins. In order for the spin diffusion process to take place, the nuclear 
quadrupole energy difference must be taken up by another thermodynamic reservoir, for 
example, by the dipole-dipole energy one. 

Recentl y th e theories for spin diffusion of the nuclear dipolar order via Pi's 0,0] in NMR 
and NQR and spin diffusion in an inhomogenius feild have been developed. It was 
shown that thermodynamics resevoir of the dipolar order plays an importent role in the 
spin latice relaxation and spin diffusion in an inhomogenius feild. Nuclear dipolar order is 
characterized by a state with nuclear spins oriented along an internal local field generated 
jy the dipole-dipole interactions (DDI) and can be described by a dipolar temperature 
1^ 12, 14, Here we consider the phenomena of spin lattice relaxation and spin diffusion 
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for both the nuclear quadrupole and dipolar energies of the nuclear spins due to their DDI 
in solids containing Pi's in an inhomogeneous field. 
Theory 

1. Hamiltonian 

The evolution of the spin system consisting of nuclear spins with / > 1/2 and PI spins 
may be described by a solution of the equation for density matrix p (t) (in units of h = 1) 

^^ = [nit),p{t)] (1) 

with the Hamiltonian 

n{t) =nQ + Hm + Hpi + np + Ubr (t) . (2) 

Here Tig represents the interaction of the J-spin system with the EFG; Tidd and Tipj are 
the Hamiltonians of the dipole-dipole interaction between nuclear spins and nuclear and PI 
spins, respectively; Hp describes the impurity spin system; Hbr (t) = Yll=-2 E^~'^'' (t) the 
spin-lattice interaction Hamiltonian , describes spin-lattice relaxation caused by the torsional 

n 

vibrations (Bayer mechanism) |lq , where A"^ is a bilinear function of the spin operators and 
(t) is a random function of tim e |l4l| . 
Using the projection operators |l7| e(^„ and defined by their matrix elements 
(m' \ei^n\n") = S^'m^n'n ^^"i W kycrlcr') = 6j^'^,6^'^ and introducing a projection density 
operators, Cmn (r), for the nuclear spins /, and Emn (r) for PI spins 

emn (r) =^Hr- r^) e^;„; (r) = ^ 5 (r - rj) ef„„ (3) 
M j 
the density of the Hamiltonians Tig, Tidd , and Tipj can be written down in the following 

form: 

Hq (f) = (2/ +l)-'Yl ^'mnemm (r) , (4) 
mn 

(r) = / dr- sin (r- ^) e„„ (r) e„v' {^) , (5) 

mnm n 

-Hpi (r-) = [ Yl ' (^"'- (^1 ^m'n' , (6) 



g mn 

where = Am — A„, Am, |^) , and \n) are the eigenvalues and eigenvectors of the operator 
Hq. {u) and \a) are eigenvectors of the operator Hp; matrix elements gl^^ {r — r') and 
/mn" ^) can be presented as 



with ft = n and 



(r - r-^) = Fmn" (r - r-') {6mn + 6,,) {6mn + S,,) , (9) 

/ / / / 

cind -Fmn" ^^^^ matrix elements of the dipole-dipole Hamiltonians Hdd and Hpi in 



Tin -represent at ion 



id. 



2. Diffusion equations 

To obtain the equation describing the spin diffusion and spin lattice relaxation of both 
the quadrupole and dipolar orders we will use the method of nonequilibrium state operator 
[isj ]. which has been applied to obtai_n the diffusion equation in cases of the Zeeman order 



spin diffusion [l^ and dipolar order ^ spin diffusion. 

Using the commutation rules between the components of the projection operators (3) 

[Cmn (r) , am'n' (^)] = 5 (f- r") (^n^'Cmn' (^) - Sn'mCm'n (r)) and [Cmn (r) , (f )] = 0, We 

can obtain the following equations in the form of localized laws of conservation of the spin 
energy densities 



dt 



+ div [jmm {r)j = Kmm (r) + Lmm (r) , (10) 



+ dtv yjdd {r}) + 2^ Jmm [r) — = KdP (r) + LdP (r) , (11) 



at \ / or 



dt / \^ ' dt 



^ ^ Pmn (^) ^mn (r) ^mm {r) 



+ — } ■ (12) 



The last equation is result of the energy conservation law. In Eqs (10) jmm (^) is the flux of 
operator e^m {r), 



4 



r = -- 



efcm (^) e^fc (^) + efcm (r) e^^ (f )} 



(13) 



-ft'mm (^) = «[^p/, Cmm {f)\ and L^^ (f) = i\Khr, {t)] in Eq. (10) are the change of the 
nuclear quadrupolar energy density due to the interaction with the PI and caused by the 
torsional vibrations (Bayer mechanism) , respectively. In Eq.(ll), jdd {r) is the operator of 
the flux of nuclear dipolar energy, 



Jdd (r) = -2z ^ j j {^ " " ^) {^ " " (^') (^) 



mnpq rn n -p 

Gil I ( r — rie ' (f) — G i i i ( r — r]ei (r) 

^ m n qp \ I ^pp \ J ^ m n p p \ I p q \ J 



(14) 



Kdp (r) = i [Hpi, Hdd {r}] and LdP (r) = i [Hbr, 'Hdd {r)] in Eq. (11) are the change of the 
nuclear dipolar energy density due to the interaction with the PI and caused by the torsional 
vibrations (Bayer mechanism) , respectively. Note, that in the case with a homogenous 
magnetic field, ^"^^^^^ , from the system equations (10) and (11) we have two separate 
equations: Eq. (10) leads to the localized law of conservation of the quadrupolar energy 

n n 

densities |20[ and Eq. (11) leads to conservation law of the dipolar energy [10[. 

In the high temperature approximation we can write the density matrix in the following 
form Q 



cr 



- ^ dX[B{t + i\) -{B{t + zA))] I aeq, (15) 

where the thermodynamic average (...) corresponds to an average with the quasi-equilibrium 
operator aeq=e~-^/Tre~-^, and 



A 



df 



{2I + 1)-'J2 Z^-" <n ^rnm (r^ + Pd (r) Hdd (r) 



PpHi 



(16) 
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B{t + tX) = e-^^ I dte'' I df{21 + 1) ' JZ^^-- ^) '^^-^ ^) + ^) ^ 

J —OO J ™„ 



X (f, t) - (3, (f, t)] ^^^^^^ + J, (f, t) VP, (r, t) 



+ (r, t) - i^zs (r, t) + [/S^ (f, t) - Pl] K^s (r, t)}e^^. (17) 

By using Eqs. (10) -(16), and taking into account that for single crystal sample of cu- 
bic symmetry, the diffusion coefficients, both for quadrupolar and dipolar energies, which 
in the general case of noncubic symmetry is a symmetrical tensor of second rank B], re- 
duces to a scalar quantity. By introduction the quantities [[3mn {r, t) — Pl] = S,mn if^ t) and 
[Pd {r, t) — Pl] = C (^; t)y fhe diffusion equations can be obtained 

d^mn {r, t) 1 



dt uj9„„ (r 



- V {D^n (r) K„ (r-) VUn (r, t) + {Un (r, t) - ( (f, t)) V^l (f)] } - 



mn 



-Wmnif^ar,t), (18) 

^ ^ [u;'_ (rl VUn (r, t) + {Un (r, t) - ( (F, t)) (r-)] + 

(1/6 IVlfjifYi 
mn 

+ V [D, (f) VC (r, t)] - W^, (f) C (r, t) . (19) 

where M^m = J dfC^ (f - r") (e^„ (f)). 

The first term in the square brackets of the right side of Eq. (18) describes the time 
dependence of the (^5 1) due to the spin diffusion with a diffusion coefficient of 

Dmn (r) = l^xj dte'' (Jmn (r. A, t) jmn (r")) /Tr {el^ (f)) (20) 

The second term gives the variation of Pmn (^, t) due to interaction with the dipolar reservoir 
in the inhomogeneous field. The last term in the right side of Eq. (18) gives the relaxation 
of (^; t) toward the inverse lattice temperature with density of the transition probability 
per unit time, Wmn (^) , which for a cubic crystal is given by {Kmm (^) + -^mm i^)), where 

/•l /"OO 

Wmn {r) = dX dte"^ {{K^m (r. A, t) + L^m (r. A, t)) {K^m {r) + L^m (r))) /Tr {emmCmm (r)) 



J-OD 

(21) 
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The first term in the curly brackets of right side of Eq. (19) describes the time variation of 
the dipolar energy due to the spin diffusion with diffusion coefficient 



D, (r) 



1 /-oo 

d\ 

^ -oo 



dte'' U (r, A, t) J, (f) /Tr {HdHd (F)) 



(22) 



The second term gives the variation of (3d (^, t) due to the interaction with the quadrupole 
reservoir in an inhomogeneous field. The last term in the right side of Eq. (19) gives the 
relaxation with the density of the transition probability per unit time, W^, (r), which for a 
cubic crystal is given by 

/•I /•oo 

Wa (r-) = d\ dte'' {[K^p (F, A, t) + L^p (F, A, t)] [K^p (f) + L^p (f)]) /Tr {HMd (r)) 



(23) 

The boundary conditions can be introduced by defining a sphere with radius I about each 
PI, called the spin diffusion barrier radius. Inside this sphere the spin diffusion process goes 
to zero: 

y(3mn (r, t) I \^=i= and V (3^ (r, t) \ \^=i= 0. (24) 

To obtain the radius of spin diffusion barrier, let us emphasize that -Dmn (0 and (r) are 
a function of the distance r from the nearest PI. In the Gaussian limiting case the stochastic 
theo.. ot .a.ne«c R d.— dependence of the d.ffn.on coeffi.en. 

Dmn if) and Dd if) can be expressed by the following function [3] 



{r) , Dd (r) 



exp 



V<„ (r) f 



Using Eq. (26) the diffusion barrier radius 



111, Q, y, y , Q 



(25) 

for the spin diffusion of 



quadrupole energy can be found by solving the equation. 



37s 

/3 



ge/i 



i± 1 



677 



|3 ' 



(26) 



where the first term in square brackets in Eq. (26) describes a distortions of the crystal 
field as a result of the inclusion of the PI. In Eg. ( 26) ro is distance between neighboring 



2J|. It was assumed that the distortion of 



nuclei, ^ is the Sternheimer antishielding factor 
the electric field is equivalent to the presence of a charge q 

The examination of the functional dependence (25) for the dipolar diffusion coefficient 
results that the main term for Dd (r) does not include any dimensional dependence. Thus 
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the radius of the diffusion barrier for dipolar energy is = |ro| , which corresponds to 
non-barrier diffusion and to the fastest relaxation of the dipolar energy. 

In the case of a homogeneous magnetic field, Vcu^^ (r) = 0, equations (20) and (21) 
give the results obtain earlier for the spin diffusion of the quadrupole 5] and of the dipolar 
energies |8i]. From Eqs. (20) and (21) we get that the dissipation of the density quadrupolar 
and dipolar energies are driven by: i) the exchange between them; ii) spin diffusion process: 
and iii) direct relaxation to the PI. 

2. Direct relaxation regime 

Exact solutions of Eqs. (20) and (21) are extremely difficult, even for simple a model 

situations. That is why we consider evolution of the spin system in time by using the next 

considerations. Immediately after a disturbance of the nuclear spin system, the gradients 

of ^mn t) and ( (r , t) are sufficiently small and diffusion cannot be of importance at the 

I I 

start of the relaxation process '22], this is the so called diffusion vanishing regime [? ]. To 
describe the relaxation at that time interval we can use Eqs. (20) and (21) by putting all 
inverse temperature gradient-terms equal zero, VC,mn{^,t) = and VC(^, ^) = 0. We also 
accept the approximation that at distances larger then the radius of the diffusion barrier 
the diffusion coefficient is independent of r 0]. Under these approximations Eqs. (20) and 
(21) come to 

d^mn {r, t) D 



dt (r) 



{[{Un (r, t) - C (r, t))]} - Wmn C [t, t) , (27) 



'■mm 



The evolutions of the C,mn {"f^yt) and C{r,t) toward their steady-state values is a linear com- 
bination of (2/ + 1) exponents with = (2/ + 1) relaxation times tn (r)- These relaxation 
times r^v {r) are the function of the position r. In order to obtain the experimentally ob- 
served signals, the solutions of Eq. (27) and (28), must be averaged over the sample. For 
this averaging procedure one needs the knowledge the field distribution. 

As a result of the diffusion vanishing relaxation regime the local inverse temperatures, 
^mn{'^,t) and ({r,t), become spatially distributed over the sample with a distribution which 
is not the equilibrium one. In this case we have to take into account also the gradient-terms 
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, ^Cmn{'i^,t) and V({f,t) in Eqs.(20) and (21). In the next section we will consider the 
influence of the spin diffusion process. 

3. Diffusion relaxation regime 

Assuming that at distances larger then the radius of the diffusion barrier the diffusion 
coefficient is independent of r ^j. Multiplying Eq. (20) by 'yoo'^.^ (r) (cmm {r)) and Eq.(21) 
by J2m and by integrating Eqs. (29) and (30) over space variable r , we obtain 

the equations describe the evolution of the experimentally observed values, the quantaty 
that conects with the Z component of the total nuclear magnetization [Mz (t) — (0) = 

Xlmn -^mn (t)) ^ -E'mn (t) 



dEmn (t) 7 



^ /rffE,(f,t)V^L(r)] 

- J dfWmnif^Err,n{r,t), (29) 



and the total dipolar energy, (t) 
dE,{f,t) 



dt ^ 



ran 



Dmn j dfEmn (r, t) Au^^ (f) - j dfWd (F) Ed (f, t) , (30) 



where Emn (t) = J dfEmn {r, t), Emn {r, t) = 7u;^„ (f) ^rnn {emm {r, t)), Ed (t) = J dfEd (r, t), 
andErf(f,t) = ^^M^^C(r,t). 

To obtain the solution of the Eqs. (29) and (30) and to calculate the relation times, 
both for the nuclear magnetization and for the dipolar energy, we need to know the internal 
distribution function of the field, uj^^{r). As it follows from Eqs. (29) and (30) that for 
a special distribution of the internal field, Au;^^ (r) = 0, the diffusion equations gives two 
uncoupled equations: 



j drWmn{^Eran{r,t), (31) 



dt 
and 

^Ml = _ I drWd (r-) Ed (r-, t) . (32) 
Solving Eqs. (31) and (32), we obtain the normalized relaxation functions 

D _ -^mn {t) - (oo) _ -t/Ti^„ (oo\ 

rimn ItJ - — — - - — — - e [66) 
J^mn (Oj — J^mn [OO) 
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and 

^ ,, Edit) — Ed (OO) ,/rp , . 

Ed (0) - Ed (oo) 

where Ti^n = — ^'1% and Tid = — ^'1% , Cp is the concentration of the PI and E is 
the angular average of the couphng dipolar constant of the DDI between the nuclear and 
PI. It follows from the solutions (33) and (34) that at a long time after the excitation 
of the spin system, the nuclear magnetization describes by the sum of exponents — 
Ylimn "^mnC"*^'^^'"" , while the dipolar energy decreases to equilibrium exponentially. 

Results and discussion 

We will cornpare the results obtained here with the relaxation processes of the nuclear 
magnetization p] and the dipolar energy in the mixed state conventional superconducting 
vanadium (J = 7/2). In the type II superconductors, an applied magnetic field Hq, in the 
range between the lower and upper critical field , Hd < Hq < Hc2, penetrates into the 
bulk sample in the form of vortices, each with a quantum flux of $o = |f , which form a 
two-dimensional structure in the plane perpendicular to Hq [7 ]. The distribution of the 
internal vortex field can be obtained by solving the Landau- Ginzburg equation, which gives 

^(p) = ^ln^forr<A (35) 

where A is the London penetration length and p is the distance from the core of vortex in 
the cylindrical coordinate , r"^ = p"^ + z"^ . Using experimental data 0, Q| we obtain the spin 
relaxation time for the nuclear magnetization, Ti = 43 sec and the spin relaxation time for 
the dipolar energy, Tid = 93 msec . As a consequence, the dipolar energy decreases to the 
equilibrium state with an anomalously short time as compare to the relaxation time of the 

nuclear magnetization, ^ = 442. Theoretical estimation of the ratio of the relaxation time, 

^ 3 



using Eqs. (33) and (34), results ^ = y^^j .Taking into account that the distance between 
neihboring vanadium nuclei = 2.63 x lO'^cm and the radius / = 1.96 x 10^^ cm , for the 
ratio ^ = 413 , which is in a good agreement with the result obtained from experimental 
data. 

Conclusions 



In conclusion, we obtained coupled equations describing mutual relaxation and spin dif- 
fusion of the quadrupole energy and dipolar energy by using the method of nonequilibrium 
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state operator [18[. The equations were solved at a short and long times approximations 
corresponding to the direct and diffusion relaxation regimes. We showed that at the begin- 
ning of relaxation process, the direct relaxation regime is preferred. The relaxation regime 
changes both for the nuclear quadrupolar and the dipolar energies, to the diffusion one. 
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